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Formation of domain walls during a rapid phase 
transition in a quasi one dimensional Cahn-Hiliard 
equation describing binary fluids in a thin tube is 
studied. Density of kinks scales like a sixth root of 
quench rate for equal concentrations and like a square 
root of quench rate for unequal concentrations of flu- 
ids. For a slow inhomogeneous transition the density 
is linear in velocity of temperature front. This paper 
is flrst theoretical study of topological defects forma- 
tion in a system with conserved order parameter. 

It has been pointed out some time ago that topolog- 
ical defects which formed during subsequent symmetry 
breaking phase transitions can provide seeds for struc- 
ture formation in the early universe 0. Kibble gave 
a detailed theory of defect formation in I order phase 
transitions which proceed by bubble nucleation. Bub- 
bles are born with random orientation of order parame- 
ter; when they coalesce they can give rise to a nontrivial 
topological winding number. This prediction was veri- 
fied in relatively simple beautiful experiments in liquid 
crystals |3|. Disclinations were observed and theoretical 
relation between bubble density and disclination density 
was verified. All this at nearly room temperature and 
with at most the aid of an optical microscope. 

A scenario for II order transitions was put forward by 
Zurek jj]. He observed that the order parameter goes 
out of equilibrium some time before the critical point is 
crossed. In this way a unique finite correlation length ^ 
proportional to the forth root of the quench rate is frozen- 
in. ^ is the scale which determines density of topological 
defects after the transition. Experimental verification is 
not as clear as for I order transitions. Experiments in 
superfluid "^He were done just to be falsified later on. 
There are spectacular experiments in superfiuid ■^iJe 
Density of detected vortices is consistent with the order 
of magnitude Zurek's prediction. However the tricky de- 
tection is somewhat indirect, there is no way to control 
the quench rate and verify the scaling of ^ with the rate, 
finally the precise origin of detected vortices was also re- 
cently put in question There is urgent need for an 
experiment as simple as that in liquid crystals. 

A good candidate for such an experiment are bi- 
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nary fluids described by Cahn-Hiliard plus Navier-Stokes 
equations. Order parameter is conserved and propor- 
tional to the difference in fluid concentrations. An ex- 
ample are aniline and cycloheksan. They do not mix for 
temperatures less than Tc = 30.9°C. The transition is 
second order with a standard Ginzburg-Landau free en- 
ergy. The fluids differ in optical density and thus domain 
walls between them can be detected by optical means. A 
rapid quench should result in a quench rate dependent 
density of domain walls. The best experimental approach 
is to quench the fluids in a thin tube, which makes the 
system effectively one dimensional. In one dimension the 
quench-made domain walls cannot be eradicated by mu- 
tual annihilation thanks to order parameter conservation 
- they are a permanent record of the nonequilibrium tran- 
sition. The transition in one dimension is, strictly speak- 
ing, a crossover; the correlation length does not diverge at 
Tc- This deviation from a mean-field is important only 
in a very narrow critical regime which is irrelevant for 
our nonequlibrium quenches. In fact we do not see any 
deviation from a mean-field in our numerical simulations. 

Dynamics of binary fiuid separation in one spatial di- 
mension is described by model H dynamics which consists 
of the Cahn-Hiliard (CH) equation for conservative real 
order parameter 

^=[-e{t,x)cj, + <j>' -<!>"]" + e , (1) 

where ' = dt and ' — d^. CH equation is usually supple- 
mented by Navier-Stokes equation but in one dimension 
the incompressibility condition makes this perturbation 
trivial. We allow for variation of the symmetry breaking 
parameter e (temperature/pressure) both in space and 
in time. ^ is assumed to be a white Gaussian noise with 
nonvanishing cumulants 

{^{tl,Xi)^{t2,X2))=2T5{t^^t2)5{xi-X2) . (2) 

TRANSITION WITH EQUAL CONCEN- 
TRATIONS. To begin with let us consider a uniform 
linear phase transition with 

e{t,x) = - , (3) 

T 

with = on average at the initial t = —oo, which is 
preserved by CH evolution. Any uniform transition close 
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to the critical point at e = can be described by such a 
hnearized e. 

For t < the system is in a symmetric phase; (f) is 
subject to smah fluctuations around 0. For this stage of 
the quench and also for the onset of spinodal instability, 
just after e crossed 0, Eq.(|l|) can be linearized in 0, 
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The mean-field (0^) is kept just to control validity of the 
linearization. We solve Eq.(^) by neglecting the mean- 
field term and performing a Fourier transformation 
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Eq.(Q) can be solved for any k with a help of its Green 
function and the correlations (^. The power spectrum 
P{t, k) of the fluctuations is 



b*{t,k)4>{t,p)) = p{t,k)S{k-p) 
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The fluctuations, measured by 
small for t < 0. At some t > 



J dk P{t,k), are 
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and the linearized approximation involved in solving 
Eq.(^ breaks down. Assuming that tjT^I'^ > is suffi- 
ciently large the error function integral in is constant 
for small k^ (but it is still essential to suppress the diver- 
gence for large fc^). The remaining exponent is peaked at 
k^ K, t/Zr with a maximum of exp(4t^/27T^). This maxi- 
mum begins to blow up at £ « t"^/^ . It is around this time 
that {(fP') passes through t/r. This is the moment when 
kinks of width given by the corresponding healing length 
of T^/^ begin to form. Fluctuations with |fc| > fc = r^^/® 
are irrelevant for kink formation at t. Density of kinks 
which are going to form can be identified with the den- 
sity of zeros of (j){t,x) smoothed over |fc| > k which is, 
according to a formula from 



of adiabatic transition r — > cx) the system stays close 
to the critical point for a time long enough to order at 
long distance. For a fast quench substantial amount of 
initial disorder is frozen into the ordered phase in a form 
of kinks. In one dimension, because of ^-conservation, 
these domain walls have no freedom to diffuse around 
and mutually annihilate. There is no "phase ordering 
kinetics" to erase this trace of disorder. The kinks are 
permanent record of the phase transition. 

TRANSITION WITH UNEQUAL CONCEN- 
TRATIONS. The two fluids may differ by average 
concentration. In that case the conserved average </> is 
M 0. We take M > for definiteness. An uniform 
(f) — M configuration is stable against small perturba- 
tions if e < 3M^ when M is outside the interval between 
the two inflection points of the double-well potential. For 
e < (j) — M IS, not a minimum of the effective po- 
tential but it cannot decay because of ^-conservation; 
<f> = M is bigger than the positive minimum at ^/e for 
e > or the minimum at for e < 0. 

At e = (j) = M coincides with the minimum 

of the potential. When AP < e < 3A'P = M is 
stable against small perturbations but again it is not a 
minimum of the double well potential. This time, how- 
ever, 0-conservation does not forbid its decay to the min- 
ima at ±-\/e. The decay can proceed by thermal nucle- 
ation of antikink-kink pairs (AKP) provided that T is 
large as compared to a barrier. To estimate minimal T 
(7) necessary for nucleation let us introduce e = + e 

with e < and expand = M + (j>. The effec- 
tive potential can be approximated by i|e|0^. Fluctu- 
ations around M are {(jP) ~ T/ ^/\s\. They can re- 
sult in AKP nucleation if they reach beyond the posi- 
tive inflection point at -y/ (SM^ — |e|)/3 or in other words 
{M - ^(3Af2 - |e|)/3)2 « {(f>^). To first order in \e\/AP 
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Introducing an integration variable k/k one can see that 
any r-dependence can be factorized in front of the inte- 
grals so that 
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for any r. Results from numerical simulations consistent 
with this prediction are shown in Fig.l. In the limit 



nucleation takes place at |e| < {TM ) '° = e„. Nucle- 
ation time is ~ e"^. If the transition proceeds at a finite 
rate, e = t/r, there may be not enough time for thermal 
nucleation. Nucleation can only happen if at e = — e„ the 
nucleation time is shorter than the time left till e reaches 
0. This condition is satisfied if 
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If the transition is fast enough or T is sufficiently small, 
when the opposite condition holds, no AKP are nucleated 
for IvP < e < SM^; (p remains fluctuating around M 
until e crosses 3M^ and spinodal decomposition due to 
instability oi (j) = M begins. This case can be analyzed 
as follows. As a first step we define e(t) = SM^ + t/r 
and expand 4> = M + cf). Equation (|l|) when linearized 
in (f) gives Eq.(P but with (f) replaced by (f>. By the same 
token as for equal concentrations, at f = r^/^ fluctuations 
blow up with momentum peaked at ±A; ~ ±r~^/^. In 
contrast to the M = case, the growth of this instability 
is halted very quickly much before any kinks are formed. 
To see this let us pick a wave (j) = o,(t) cos{kx). At around 
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i = T^l'^ maxima of this wave enter the area beyond 
the inflection point of the actual double-well potential, 
> ^ (3Af2 + T~^/3)/3 - their growth is slowed down. 
At the same time minima get further and further into 
the unstable regime between the inflection points so it 
seems that their growth rate should be accelerated. If 
it were so they would quickly hop to the neighborhood 
of the negative inflection point and the initial cos(fca;) 
would distort into a periodic array of kinks-antikinks with 
density ~ fc. It cannot be so because it would obviously 
violate ^-conservation. The growth of a(t) stops as soon 
as the maxima cross the inflection point. At this stage 
the cosine wave is still an almost negligible fluctuation 
around the large M, especially for large r. However the k 
scale is frozen-in as a record of the quench. That is what 
remains of the spinodal stage. AKP nucleation begins. 

The problem with the cosikx) lies with its perfect pe- 
riodicity which cannot be broken. However the power 
spectrum although peaked at ~ ±fc is also spread around 
them by ~ fc. This means in particular that some min- 
ima of the "cosine" are more negative than the other. 
There is competition between them. Each of them wants 
to hop to the negative inflection point but the number 
of winners is limited by ^-conservation - the average 
must remain M . 

A fraction of winners can be estimated in a follow- 
ing way. At the end of the spinodal stage we have a 
landscape of correlated domains of size hr^ . A frac- 
tion q of negative domains will win. Those negative do- 
mains who win will end at the negative inflection point 
(f> = — \/ + 1/3t^/^. All the positive domains and the 
negative loosers will end at = +\/M^ + 1/3t^/^. To 
keep 4> — M on average, q must satisfy 

M = ['^+(^^'^)] + ^ yAf^ + l/3rV3 , (11) 

where [—q+{l — q)] comes from positive and 1 from nega- 
tive domains. For M^r^/^ > 1 we obtain q w 1/6]VPt^/^ 
so the kink density should scale as 

" - 'Jfc-' - for » 1 . (12) 

This 1/2 exponent is 3 times bigger than the 1/6 in 
Eq.(^. It should be much easier measurable in exper- 
iment. Numerical simulations are consistent with this 
prediction as shown in Fig.l. The process of spinodal 
decomposition followed by nucleation is shown in Fig. 2. 

INHOMOGENEOUS TRANSITION. An uni- 
form phase transition may be a good flrst approxima- 
tion in some cases but in real life we often have to face 
the fact that it is not perfectly homogeneous. To gain 
some insight let us begin with the temperature front of 
the form 

N / +e- , X <vt , . 

^(^'^)^\-.+ , vt<x (13) 



in the limit of very slow v. We will argue a posteriori that 
the sharp step is a good approximation of any generic 
front for w — > 0. 

We solve the problem of kink generation behind the 
moving front by perturbative expansion around v = 0. 
At w = there is a static 0-front, (j){x) = H{x) - a step 
in at a; « interpolating between — y/eT aX x ^ ~oo 

and at a; = -|-oo. Its width is given by w e_ + e, . 

Let us now switch on small w > 0. The e-front (O) is 
slowly moving on. If were not conserved, the 0-front 
would follow moving in step with e-front and leaving no 
kinks behind For our conserved 4> it is not possible, 
kinks must inevitably appear. To see in some detail how 
it happens let us substitute (j){t, x) = H{x — vt) + ip(t, x) 
with — 0{v) to Eq.(0) and keep only 0{v) terms. We 
are interested in length scales large as compared to the 
width of the step H{x) and that of the e-front. That is 
why we keep only up to the second x-derivative which is 
responsible for diffusion. Far from x k vt we obtain 

^{t, y) = e+r{t, y) + v^'{t, y) + e{t)5{y) , (14) 

where y — x ~ vt. We take into account a i5-like source 
term at y = which is a long-wavelength approximation 
to vH' . We also set e_ = e+/2 for simplicity. The source 
term is switched on at t — when the e-front starts 
to move, hence the Heaviside function. The solution of 
Eq.(|l4[) is straightforward, 

[„+„(t_t')]2 

m y) = « V f / dt' . (15) 

The source term produces ip ai y = Q aX constant rate. 
It spreads around by diffusion but at the same time it 
is carried to negative y with velocity —v. For y > 
the penetration by diffusion dominates at first but at 
v'^t^ ^ e^t the two processes balance one another and 
ip{t,y > 0) saturates. From this time on all the ip is 
carried directly from the source to y < with velocity 
—V. This effectively means that the ip-iront halted, while 
the e-front keeps moving on. A supercooled phase with 
a slightly positive (j) is growing in between them with 
velocity v. When its width exceeds A/2/e+ , (p decays 
towards positive ground state. From this time on we 
have a negative 0-step moving together with e-step and 
the whole story repeats itself at spatial intervals of e_(_/w. 
Density of kinks is 

v 

n ^ — for u — > . (16) 

It should be stressed that the whole process is determinis- 
tic, kinks are made at regular intervals. Noise is required 
to begin the process; it also adds some irregularity on top 
of the regular pattern. 

Note that for small v the relevant length scale is e+/v. 
For small enough v, it far exceeds the e-front width and 
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the width of H{x). This justifies the sharp step in Eq.jl^) 
and the long-wavelength approximations involved in our 
derivation of Eq. (|lj) . 

Let us now turn to the opposite large-u limit where we 
anticipate the transition to be effectively homogeneous. 
Any generic e(t,x) can be linearized around e = 0, 

vt ~ X 

e(t,x) — = a{vt — xj . (17) 

VT 

At any fixed x the transition proceeds at the rate of 1 /r 
just like in Eq.(|3|). If it were homogeneous it would en- 
hance the momentum k — r^^/^. For 

> or w > r"^/^ (18) 

this momentum scale is much bigger than the slope a 
and the relevant field fluctuations do not feel the inho- 
mogeneity. This is where the transition is effectively ho- 
mogeneous and Eq.(^ applies. 

CONCLUSION. We studied dynamics of domain 
walls formation during a quench in effectively one dimen- 
sional binary fluids. We believe that our predictions can 
be verifled by experiments in thin test tubes. In fact first 
steps in this direction has been already made jl^. One 
can test the 1/6 scaling for equal concentrations, which 
is the most direct analogue of Kibble- Zurek scenario for 
second order transitions. A more intriguing case are un- 
equal concentrations where we also get scaling with the 
quench rate in a range of parameters. Its 1/2 exponent 
should be much easier to measure than the 1/6. Finally, 
unlike for nonconserved order parameter Q, density of 
kinks is linear in front velocity for a slow inhomogeneous 
quench. We believe that binary fiuids are a unique op- 
portunity of an almost bare eye detection of topological 
defects. We also believe that one can not only detect 
topological defects but also proceed and measure scal- 
ing of their density with quench rate. In this paper we 
lay theoretical foundations for the experimental work in 
progress [ p^ . 
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FIG. 1. log(n) as a function of log(r) for M = (the top 
plot) and for M = 2 (the bottom plot) according to numerical 
simulations. For A'l — the slope is 0.18 ± 0.01 as compared 
to the theoretical 1/6 ~ 0.17. For M = 2 the slope saturates 
for log(T) > 2 at 0.48 ± 0.05 as compared to the theoretical 
0.50. At low T the M = 0, 2 results tend to be the same. 
Vertical size of a point is its statistical error. Simulations 
were done at T = 10"^ on a Ax = 1, At = 0.01 lattice of 
size 1024 with periodic boundary conditions. e(t) was swept 
from 3M^ - lOr"^/^ to 3M^ -I- lOr"^''^. Kinks were counted 
at final time. Density n is an average over many runs. 
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FIG. 2. Two snapshots of </> as a function of x for 
M = 2,T = 128, T = 10~® taken from numerical simulations. 
The thin line is </> at t = 0.8t when spinodal decomposition 

begins. The amplitude of fluctuations around M = 2 is mag- 
nified 100 times. The thick line is at i = l.Qt when kinks 
are already well defined. 



